function Z=jianmo(t,F)
V=F(1);
theta=F(2);
x=F(3);
z=F(4);
m=F(5);

S=pi*3.35*3.35/4;
g=9.8;
Isp=2958+30*randn(1);
ma=V/340;
Ma = [0.6 0.8 0.9 0.95 1.05 1.1 1.2 1.5 2 2.5 3 3.5 4 4.5 5];
H1=0:20000:80000;
zhi=[2.926	3.237	3.445	3.571	4.111	4.614	4.581	4.815	4.800	4.723	4.673	4.623	4.571	4.520	4.469
    3.080	3.407	3.626	3.759	4.327	4.857	4.822	5.069	5.053	4.972	4.919	4.866	4.812	4.758	4.704
    3.234	3.577	3.807	3.947	4.544	5.100	5.063	5.322	5.305	5.220	5.165	5.109	5.052	4.996	4.939
    3.696	4.088	4.351	4.511	5.193	5.828	5.786	6.083	6.063	5.966	5.903	5.839	5.774	5.709	5.645
    4.466	4.940	5.258	5.451	6.275	7.042	6.992	7.350	7.326	7.209	7.132	7.055	6.977	6.899	6.821];
if ma<0.6
    ma=0.6;
elseif ma>5
    ma=5;
end
Ca=interp2(Ma,H1,zhi,ma,z);
Ca=Ca*(1+0*unifrnd(-0.1,0.1,1));
Ma=[0.6 0.8 0.9 0.95 1.05 1.1 1.2 1.5 2 2.5 3 3.5 4 4.5 5];
cn=[1.966 1.902 2.143 2.244 3.3 3.624 3.33 3.307 3.256 3.021 2.652 2.365 2.154 2.026 1.897];
if ma<0.6
    ma=0.6;
elseif ma>5
         ma=5;
end
Cn=interp1(Ma,cn,ma);
Cn=Cn*(1+0*unifrnd(-0.1,0.1,1));
load('U.mat');

    P=interp1(U(3,:),U(1,:),t);
    alpha = interp1(U(3,:),U(2,:),t);

rou0=1.225;
beta=0.0001208;
Z=zeros(5,1);
rou=rou0*exp(-beta*z)*(1+0*unifrnd(-0.05,0.05,1)); %%%%%
D=0.5*rou*(V+0*randn(1))^2*S*Ca;
Y=0.5*rou*(V+0*randn(1))^2*S*Cn*alpha;
Z(1)=-P*cos(alpha)/m-D/m-g*sin(theta);
Z(2)=-P*sin(alpha)/m/V+Y/m/V-g*cos(theta)/V;
Z(3)=V*cos(theta);
Z(4)=V*sin(theta);
Z(5)=-P/Isp*(1+0*unifrnd(-0.03,0.03,1));
end

    

